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■^ ' We investigate the dependence of thermodynamic properties of black holes on the choice of 

^^ ■ statistical ensemble for a particular class of Einstein-Maxwell-Gauss-Bonnet black holes with cos- 

Tij- i mological constant. We use partial Legendre transformations in the thermodynamic limit in order 
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en '_ to compare the results in different ensembles, and show that the phase transition structure depend 

on the choice of thermodynamic potential. This result implies that thermodynamic metrics which 

p\ , are partially Legendre invariant cannot be employed to describe black hole thermodynamics, and 
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partly explains why a particular thermodynamic metric has been used so far in the framework of 

black hole geometrothermodynamics. 
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I. INTRODUCTION 



In the Euclidean path-integral approach to black hole thermodynamics, originally de- 
veloped by Hawking et al. [11-14], the thermodynamic partition function is computed from 
the path integral in the saddle-point approximation, obtaining as a result the laws of black 
hole thermodynamics. Originally, only the microcanonical ensemble was investigated [l|, ^ 
because the canonical ensemble led to difficulties related to the stability of the black hole 
under consideration. Later on, the canonical ensemble was investigated by York et al. pMlOf 
by using appropriate boundary conditions. It turns out that the Euclidean path-integral ap- 
proach is not well defined until the boundary conditions of the system are properly defined. 
The various ensembles available to the system are defined by identifying the corresponding 
boundary conditions. York's approach was used to analyze other ensembles ll| and to study 
particular systems like the charged black hole in the grand canonical ensemble 10|, black 



holes in asymptotically anti-de Sitter spacetimes 12|, ll3|, and black holes in two and three 
dimensions [12|, |lj| . 

Under certain circumstances, the results obtained by using the path-integral approach 
turned out to depend on the boundary conditions 12, ll5|, [l6|. The stability of black holes 
also turned out to depend on the choice of boundary conditions and, consequently, on the 
ensemble 17|]. In fact, it was found that in one ensemble the black hole can never be stable, 
independently of the values of the black hole parameters, and in another ensemble the black 
hole is almost always stable. This result shows that the stability properties of a black hole 
are drastically infiuenced by the boundary conditions that determine ensemble. Moreover, 
it follows that this behavior must be also present in the thermodynamic limit in which 
different ensembles correspond to different thermodynamic potentials, related by Legendre 
transformations. On the other hand, an important characteristic of a black hole is its phase 
transition structure which is closely related to the stability properties of the system. It then 
follows that the phase transition structure can, in principle, also be ensemble dependent. In 
the thermodynamic limit, this would imply a dependence on the thermodynamic potential. 

The formalism of geometrothermodynamics (GTD) has been proposed in 18| as a Leg- 
endre invariant approach to describe thermodynamics in terms of geometric concepts. One 
of the conjectures of GTD is that curvature singularities of the equilibrium space are related 
to phase transitions of the system. One can then wonder whether GTD is able to handle 



the dependence of phase transition structure on the statistical ensemble. This is the main 
goal of the present work. For the sake of concreteness, we will consider here a particular 
black hole configuration in the Einstein-Maxwell-Gauss-Bonnet (EMGB) theory. We will see 
that in fact the phase transition structure, as dictated by the behavior of the specific heats, 
drastically depends on the thermodynamic potentials which are related by partial Legendre 
transformations. Consequently, the metrics that in GTD are invariant under partial Leg- 
endre transformations cannot be used to describe in an invariant manner the properties of 
such black hole configurations. 

This paper is organized as follows. In Sec. HIl we present a brief introduction into 
the formalism of GTD. In Sec. IIIH we review the main aspects of a particular spherically 
symmetric black hole in the EMGB theory with cosmological constant. In Sec. IIV[ we 
apply the formalism of GTD with a particular thermodynamic metric to reproduce the 
thermodynamic properties of the black hole configuration. Finally, in Sec. |Vl we discuss 
our results. 

II. GEOMETROTHERMODYNAMICS 

In equilibrium thermodynamics, to describe a system with n thermodynamic degrees 
of freedom, one needs a thermodynamic potential $, a set of n extensive variables {-E"} 
(a = 1, ..., n), and the corresponding dual intensive variables {/"}. Classical thermodynamics 
is invariant with respect to a change of thermodynamic potential $ — t- $ which is defined 



by meas of the Legendre transformations 19 1 



{z^}^{z^} = {^^E\n, (1) 

$ = <! - 5kiE^P , E' = -f , E^ = E^, r = E\ P = P , (2) 

where iU j is any disjoint decomposition of the set of indices {1, ...,ri}, and k,l = 1, ...,i. 
In particular, for z = we obtain the identity transformation. Moreover, for i = {1, ...,n}, 
Eq.(l2]) defines a total Legendre transformation, i.e., 

$ = $ - dabE""? , E" = -r, P = E°- . (3) 

Notice that in order to apply a Legendre transformation to a tensorial object in a particular 
coordinate system, it is necessary to use the corresponding matrix representation which can 



be computed by considering al 
a {2n + 1)— dimensional space 



the coordinates as independent, that is, as coordinates of 
18|. It then turns out that to investigate the mathematical 
structure of thermodynamics in general it is necessary to use contact geometry which is 
constructed as follows. Consider a (2n + 1)— dimensional differential manifold T together 
with its tangent manifold T(T). Let V C T(T) be an arbitrary field of hyperplanes on T. 
It can be shown that there exists a non-vanishing differential 1-form on the cotangent 
manifold T*{T) such that V = kerO. If the Frobenius integrability condition O A dQ = 
is satisfied, the hyperplane field V is said to be completely integrable. On the contrary, if 
A dQ 7^ 0, then V is non-integrable. In the limiting case A ((i0)" ^ 0, the hyperplane 
field V becomes maximally non-integrable and is said to define a contact structure on T. 



The pair (T, 0) determines a contact manifold 20|. Consider G as a non-degenerate metric 
on T. The set (T, 0, G) defines a Riemannian contact manifold. 

Let us choose the coordinates of T as Z^ = {$, E'^, J"} with A = 0, 1, ..., 2n. According 
to Darboux theorem, the contact 1-form can be written as 

= rf<l> - 6abrdE' , 6ab = diag(l, 1, ..., 1) , (4) 

where we assume the convention of summation over repeated indices. Under a Legendre 
transformation, the contact 1-form transforms as — ?■ = (i<l> — 5abl'^d&. Consequently, 
the contact manifold (T, 0) is a Legendre invariant structure. If we now impose Legendre 
invariance on the metric G, the Riemannian contact manifold (T, 0, G) is Legendre invari- 
ant. Any Riemannian contact manifold (T, 0, G) which satisfies the condition of Legendre 
invariance is called a thermodynamic phase manifold and constitutes the starting point for 
a description of thermodynamic systems in terms of geometric concepts. 

Notice that to construct a concrete phase manifold we only need to specify the metric G. 
It turns out that Legendre invariance does not fix completely the metric G. All the metrics 
we have found so far can be classified as invariant under total Legendre transformations 

G^ni = (d$ - IJE'^f + A {UE'^I') [xcddE^dl') , (5) 

or invariant also under partial Legendre transformations 

G''' = (rf$ - IJE'^f + A [EJaf^+ME'^dr , Ea = 6abE' , h = Sabl' . (6) 

Here A is an arbitrary Legendre invariant real function of Z^, and ^at and Xab are diagonal 
constant matrices that can be expressed in terms of the Euclidean and pseudo-Euclidean 



metrics Sab = diag(l, ..., 1) and rjab = diag(— 1, 1, ..., 1), respectively. The choice ^ab = Sab, 
Xab = Sab leads to the metric G^ which has been used to describe systems with first order 



phase transitions [18|, l2l|. The alternative choice ^ab = Sab Xab = Vab, which is denoted 
by G^^ , turned out to describe correctly second order phase transitions especially in black 
hole thermodynamics 2ll423l|. Moreover, in the metric G^^ the additional choice C,ab = 



I {Sab — Vab) can be used to handle also the thermodynamic limit T — t- 0. 

In classical thermodynamics, a particular thermodynamic system is completely deter- 
mined by its fundamental equation $ = ^{E"") which determines an n— dimensional surface 
in T. In GTD, we use the same idea to construct the manifold which should describe a 
particular thermodynamic system. However, some technical details must be considered in 
order to preserve the Legendre invariance of T. To this end, consider a (smooth) embed- 
ding map if : £ ^ Ti where £ is an ra— dimensional submanifold of the phase manifold 
{T,Q,G). If we consider the set {E""} as the coordinates of £, then the embedding map 
reads if : {E"} i — > {Z^(E")} = {$(E"), E^ J"(E")}. In this manner, the fundamental 
equation $ = $(£"^) appears in a natural way as the result of introducing a well-defined 
embedding map ip. Moreover, the metric G induces in £ the canonical metric g by means 
oi g = ip*{G), where ip* is the pullback of ip. The pair {£,g) is called equilibrium manifold 
if the map (p : £ ^ T satisfies the condition 



which implies that 



ip*iQ)=ip*{d^-SabPdE') = 0, (7) 



d^ = IadE\ ^ = 4. (8) 



The first of these equations corresponds to the first law of thermodynamics whereas the 



second one is usually known as the condition for thermodynamic equilibrium 24 1. 

In the case of the metric G^^ discussed above with ^ab = \ {Sab — Vab) and Xab = Vab, the 
corresponding induced metric of £ can be written as 

where Va — diag(— 1, 1, ...1). The geometric properties of the equilibrium manifold £ de- 
scribed by the metric g^^ should be related to the thermodynamic properties of the system 
described by the fundamental equation $(£'"). In GTD, it is conjectured that £ is curved for 



systems with thermodynamic interaction and that curvature singularities in £ correspond 
to phase transitions of the corresponding thermodynamic system. 

Notice that the arbitrary function A represents a Legendre invariant conformal factor. 
In practice, it does not affect the main behavior of the curvature of g^^ which is the main 
geometric quantity in GTD. Nevertheless, this function plays an important role if we demand 
that our results be invariant not only under Legendre transformations, but also with respect 



to changes of representation. In fact, one can show 25|] that a particular choice of A makes 



the metric G^ invariant under changes of representation in general. In the case of the metric 



G^^, an additional conformal factor is needed 26 1 



III. A SPHERICALLY SYMMETRIC CHARGED BLACK HOLE 

In this section, we study a particular black hole solution which is appropriate to illustrate 
the ensemble dependence of black hole thermodynamics. The Einstein-Gauss-Bonnet (EGB) 
theory is the most general theory in five dimensions that leads to second order differential 
equations, although the corresponding Lagrangian density contains quadratic powers of the 
curvature. The most general action of the EGB theory is obtained by adding the Gauss- 
Bonnet (GB) invariant and the matter Lagrangian Lmatter to the Einstein-Hilbert action 

I = K j d'x^[R + a{R^ - AR'^'^R^, + R''^^'' R^^^s) + Lmatter], (10) 

where k is related to the Newton constant, and a is the GB coupling constant. In the 
case of the Einstein-Maxwell-Gauss-Bonnet (EMGB) theory with cosmological constant, 
the matter component of the action fllOl) is given by 

Lmatter = FajiF^ — 2A , Faj3 = Ap^a — ^a,^ , (H) 

where A is the cosmological constant, and F^j^ represents the electromagnetic Faraday tensor. 
The low energy limit of certain string theories leads to the EMGB theory with cosmolog- 
ical constant; therefore, it is important to study the physical properties of exact solutions 
like black hole solutions. A particular solution was obtained in 
using the following 5D static spherically symmetric line element 

dv 

ds" = -f{r)de + T7-T + Ad^l + sin' Oi{ddl + sin^ ^s^^s)] • (12) 



27| (see also |28|-|30|) by 



The coordinate r has the dimension of length while the angular coordinates (6'i, 62) G [0, it] 
and 63 G [0, 2tt\.A 5D spherically symmetric solution in EMGB gravity with A was obtained 



by Wiltshire 3l|, using the metric ansatz f ll2p and the metric function 



^, , , r-2 r2 / 8aM 8aQ^ 4aA 

The two parameters M > and Q are identified as the mass and electric charge of the 
system. It is easy to see that the conditions for the solution (fT3|) to describe a black hole 
spacetime are firu) = and 

i + ^:^-^ + ^->0' (14) 

where th is the radius of the outermost horizon. In this work, we will limit ourselves to the 
investigation of positive a and negative definite A in order for the mass of the black hole to 
be always positive. 



IV. GEOMETROTHERMODYNAMIC ANALYSIS 

To find the fundamental equation of the black hole under consideration, we first note 
that the surface area of the event horizon is 

/ / sin^ e sin (l)ded(l)d'ilj = 27r^r% . (15) 

_ =0 i(/>=0 J 1^=0 

Moreover, the Bekenstein-Hawking en trop y is S* = -^rj- = -^ri"H^ which becomes S = rj^hj 



choosing the constants appropriately (32|. Since the black hole condition f{rH) = implies 
that 

\% - 2r|, + 2 (M - 2a) r^ - ^g^ ^ , (16) 

the corresponding thermodynamic fundamental equation in the mass representation becomes 

M = 2a + S-'^ + ^ - ^5^/3 . (17) 

Notice that to guarantee the positiveness of the mass, we must demand that a > and 

A<0. 

From the energy conservation law for black holes, dM = TdS + (j)dQ, we can derive the 

expressions for the temperature and electric potential of the black hole on the event horizon 

as 

2 3S^/-'-AS'-Q' 2Q 

^=9 s^s ' ^=-s^- (18) 



According to Davies [33|, a black hole undergoes a second order phase transition at 
those points where the specific heat C diverges. Strictly speaking, this means that we 
must introduce the concept of "heat", say Qheat, for a black hole. Of course, this is a 
problem that can be handled correctly only within the context of a physically meaningful 
statistical model which is probably the most important unsolved problem in black hole 
thermodynamics. The simplest available solution of this problem is to use the analogy 
with classical thermodynamics. Indeed, the first law of thermodynamics dM = TdS + (pdQ 
allows us to define the "heat" through the relationship dQheat = TdS so that dM = dQheat + 
(j)dQ. Then, following the standard approach of ordinary thermodynamics, we introduce the 
specific heat as 

which in this case is given by 

C -35|^ 3gl-Ag^-g^ \ 
"^ \5Q^-3Sl-ASy ■ ^ ^ 

The temperature (fT8|) is positive only in the range 3 S^^^ — AS"^ > Q^ and, therefore, the 
specific heat can take either positive (SQ^ — 35*3 — AS*^ > 0) or negative (SQ^ — 35*3 — AS*^ < 
0) values, indicating the possibility of stable and unstable states. In the limiting case 
5Q^ — 3S3 — AS"^ = 0, the black undergoes a second order phase transition during which it 
changes its state of thermodynamic stability. 

In classical thermodynamics, once the internal energy of the system is well defined, the 
analysis at the level of thermodynamic variables can be associated with a particular statis- 
tical ensemble. In the case of black holes, however, there is no unique definition of internal 
energy. For instance, the potential M{S, Q) has been associated with the microcanonical 



ensemble 34] and with the canonical ensemble |35[. In the path- integral method, in which 



the ensemble is fixed through the boundary conditions, one could also associate this poten- 



tial with the grand canonical ensemble [36|. Using this last option, we can say that the 
phase transition structure we have found above for the EMGB black hole can be associated 
with the grand canonical ensemble. Now we can consider the ensemble dependence already 
found in black hole thermodynamics. Indeed, since the stability properties of a black hole 
can drastically change from one ensemble to another, in the thermodynamic limit, in which 
the change of ensemble corresponds to a change of thermodynamic potential, the same effect 
can occur. 

8 



In the case of the EMGB black hole we are considering here, performing Legendre trans- 
formations on M[S, Q) one can derive the potentials 

H = M-(f)Q, F = M-TS, G = M-TS-(j)Q, (21) 

that are known as the enthalpy, the Helmholtz free energy, and the Gibbs free energy, 
respectively. The enthalpy 

H = M-(l)Q = 2a + S'/' - -0^^2/3 _ A54/3 _ (22) 

satisfies the first law of thermodynamics, dH = TdS — Qd(f), and can be associated with the 
canonical ensemble. Then, if we assume again that the "heat" of the black hole is defined 
by dQheat = TdS, the specific heat for fixed </> is given by 

^'^-[ &r )^~\&r)^~ 12^V3 _ 902^1/3 + 4A5 ■ ^^^) 

Taking into account the condition 125'^/'^ — O^^S^/^ — 4AS' > 0, which follows from the 
condition T > 0, we find that stable states {C^ > 0) are allowed for entropies in the range 
S"^^^ > 4m (4 — 30^). Furthermore, for 0^ > 4/3 and S > the last condition is always 
satisfied, indicating that stable states always exist in this case. Moreover, the specific 
heat fl23l) shows that the roots of the equation 125*^/3 — 9(p'^S^^^ + AAS = determine the 
locations where second order phase transitions occur. Notice that the singularities of C^ are 
different from those of Cq] consequently, the corresponding phase transition structures do 
not coincide. In the limiting case A — )■ 0, the difference is more dramatic since the specific 
heats 

indicate the existence of phase transitions in the potential M{S, Q) with no transitions at 
all in the potential H{S,(j)). This shows that the thermodynamic properties of this black 
hole are not invariant with respect to the partial Legendre transformation that relates the 
potentials M{S,Q) and H{S,(j)). 

Now we turn to the description of the above results within the framework of GTD. As 
mentioned in Sec. [TTl we must choose a metric to derive the geometric properties of the 
equilibrium manifold. In general, we have three different options, namely, G^ , G^^ and G^^^ . 
Since G^^^ is invariant with respect to partial Legendre transformations, the above result 

9 



shows that G^^^ cannot be used to describe the thermodynamics of the EMGB black hole. 
It follows that the only invariance that can be imposed in the case of black holes is with 
respect to total Legendre transformations for which we can use the metrics G^ and G^^ . In 
our experience, we have seen that G^^ correctly describes systems with second order phase 
transitions and now we want to see whether it can also correctly handle the dependence on 
the ensemble. 

In the grand canonical ensemble, the fundamental equation M = M{S, Q) is given in 
Eq.( !T7|) . Then, the coordinates of the equilibrium manifold S are E"" = {S,Q) and the 
thermodynamic potential is $ = M. The thermodynamic metric i^ can then be written as 

1 



^ = ^(3^'^'-^^'-^') 



,(3^^/=^ + AS^ - 5Q^)dS^ + dQ' 



95 
A straightforward computation results in the following scalar curvature: 

27 S^/^N{S,Q,A) 

"2" (3 54/3 - g2 _ A ^2)3 (3 ^4/3 _ 5Q2 + A s^) 
with 



(25) 



R=-T ...... ^. : J.vr::':^ .^. . ..... ^ (26) 



N{S, Q, A) = 42 Q2^^/3a - 34 SQ^A - 5 S^Q^A^ - 18 Q^S^^^ 

- 7 S^A^ + 36 5"/3A + 15 S^^^'^A^ - 162 S^ + 108 Q^^s/s _ ^27) 

From the expression for the scalar curvature it is obvious that the singularities are located 
at the points satisfying the equation 3 S^^^ — Q'^ — AS'^ —^0, which coincide with the points 
where T — )■ 0, and at the points satisfying the equation 3 S"^^^ — 5 Q^ + A S*^ — t- 0, which are 
the points where Gq — )■ oo. This proves that the curvature scalar correctly reproduces the 
thermodynamic behavior in the grand canonical ensemble. 

It would be interesting to prove explicitly the invariance of the above results with respect 
to a total Legendre transformation, i.e., when the thermodynamic potential is the Gibbs 
potential G(T, 0) = M — TS — (pQ, which is also associated with the grand canonical 



ensemble 37|. Unfortunately, it is not possible to express S and Q in terms of T and (p 
so that the Gibbs potential cannot be written explicitly. However, in the limiting case of a 
vanishing cosmological constant, it is possible to find explicitly the Gibbs potential and, as 



shown in [38( , the corresponding thermodynamic metric can be computed and the invariance 
with respect to the total Legendre transformation can be shown at the level of the scalar 
curvature. 

10 



Let us now consider the canonical ensemble with fundamental equation ( 122|) . The en- 
thalpy H is in this case the thermodynamic potential and the coordinates of the equilibrium 
manifold are E'^ = {S, 0). The corresponding thermodynamic metric is obtained from (^ as 

3 ,o A 

27" V^ 4' 



9=\^-l<P'-is'^' 



1^-2/3 (^1 _ 3^2 ^ A^2/3^ ^^2 _ ^4/3^^2 



(28) 



which leads to the curvature scalar 



R = ^(^>-^) (29) 

(125V3 - 902^1/3 _ 4A5)3(125V3 _ 9(/.25V3 + AASy ' ^ ^ 



with 



iV(^,0,A) = (12^1/3 _ 9^2^1/3 _4^^)2^ ^^4/3^7^2(4 _ 3^2^ ^2A^2/3(902_ 4^^2/3 _ 4)1 _ 

(30) 
The curvature singularities which follow from the limit 125'^/'^— 90^5'^/'^+4AS' — )■ determine 
the phase transition structure of the black hole, because they coincide with the divergences 
of the specific heat C^. The second set of singularities for which 125*^/^ — 90^5'^/'^ — 4A5' — )■ 
corresponds to the limit T — !■ and indicates the break down of the thermodynamic 
description of the black hole. This shows that the thermodynamic behavior of this black 
hole which follows from the canonical ensemble is correctly reproduced in GTD. 

V. CONCLUSIONS 

In this work, we analyzed the problem of ensemble dependence in the context of GTD. In 
the Euclidean path-integral approach it is known that the stability properties of a black hole 
can depend on the statistical ensemble. In the thermodynamic limit, a change of ensemble 
can be simply performed as a Legendre transformation that acts on the thermodynamic 
potential. It then follows that the thermodynamic properties of a black can depend on the 
choice of thermodynamic potential. 

On the other hand, the invariance with respect to Legendre transformations plays an 
essential role in GTD. So far, the metrics used in GTD can be either invariant under total 
Legendre transformations {G^ and G^^) or invariant also with respect to partial transfor- 
mations (G^^^). The important question is whether the ensemble dependence affects the 
formalism of GTD. We have shown in this work that in fact GTD can handle correctly this 
dependence and, moreover, it explains partially why in GTD a particular thermodynamic 

11 



metric must be chosen in order to correctly reproduce the thermodynamic behavior of black 
holes. In fact, by analyzing explicitly a particular black hole in the EMGB gravity theory 
with cosmological constant, we showed that the phase transition structure is not invariant 
with respect to partial Legendre transformations. This implies that the metric G^^^ cannot 
be used to describe black hole thermodynamics. 

However, the family of metrics invariant under total transformations contains two metrics, 
G^ and G^^, but only G^^ can be used for black holes with second order phase transitions. 
Practically, the difference between G^ and G^^ is only the signature: the first one is Eu- 
clidean and the second one is pseudo-Euclidean. The question is how such a particular 
distinctness can differentiate between systems with first order and second order phase tran- 
sitions. Preliminary results seem to indicate that this can be used to formulate an invariant 



definition of phase transitions 39 1 
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